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Abstract. We establish the difference between the propagation of semiclassical Wigner functions 
and classical Liouville propagation. First we re-discuss the semiclassical limit for the propagator 
of Wigner functions, which on its own leads to their classical propagation. Then, via stationary 
phase evaluation of the full integral evolution equation, using the semiclassical expressions of Wigner 
functions, we provide the correct geometrical prescription for their semiclassical propagation. This is 
determined by the classical trajectories of the tips of the chords defined by the initial semiclassical 
Wigner function and centered on their arguments, in contrast to the Liouville propagation which is 
determined by the classical trajectories of the arguments themselves. 



(N ■ 1. INTRODUCTION 

> 



The Wigner |13[ function W is a representation of the quantum density operator p as a real function 
on phase space, for a non-relativistic dynamical system with classical-quantum correspondence and 
n degrees of freedom. Here, all the underlying geometry is euclidean and we can include the Wigner 
function within the general framework of Weyl representation of operators: A —> A E C^{IR 2n ), so 
W(x) = (27ih)~ n p(x) , where x = (pi, ...p n ,q 1 , ...q n ) denotes a point in phase space. In this context, 



the quantum evolution of the density operator p t effected by an autonomous external hamiltonian 



Q_i' H can be written in terms of its Weyl [|12|] representation as 



(1) dW(x, t)/dt = {H, W t }(x) + 0(h 2 



>- . Even though the Wigner function W is generally nonpositive and thus cannot be taken as a classical 

probability density (in contrast to its lagrangian averages J W(p,q)dp , etc..) , equation (jl|) is the 
source of a somewhat widespread belief that the propagation of Wt(x) coincides in the semiclassical 
limit (h — > + ) with the classical (h = 0) propagation of a Liouville probability density. In fact, for a 
quadratic hamiltonian the semiclassical (exact) propagation of W t (x) is indeed classical. That 
is, if Xq — ► %t IS the classical hamiltonian flow generated by , then 

(2) W (x ) — > W t (x t ) = Wo(x ) . 

For a general hamiltonian H (whose Weyl symbol H does not necessarily coincide with the classical 
hamiltonian h) the classical propagation (E|) induced by (|XJ) remains a good approximation (semiclas- 
sically correct) if W is a fairly smooth function (as with the Weyl representation of an observable), 
in which case the corrections to (|l|) can be semiclassically ignored. Such a smoothness condition can 
be realized for highly mixed statistical states, in which case W will look very much like a classical 
probability density in phase space and not surprisingly its propagation will be nearly classical. 

On the other hand, pure states are not in general represented by smooth Wigner functions. Indeed, 
the semiclassical limit W for the Wigner function of a pure state ^ in one degree of freedom |Lf is: 

(3) W(x) w A^x) cos {S4x)/h - tt/4) , 
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where S^(x) is the symplectic area between an arc of a (Bohr-Sommerfeld quantized) curve ip and 
its corresponding chord centered at x , and A is a smooth amplitude function. Such an expression 
for W is highly oscillatory, specially at this semiclassical limit, and its successive derivatives knock 
off all the favorable powers of fi in the "corrections" to (|l|). The important exception is when x 
lies very close to ip in which case S^(x) is very close to zero and constant, and W(x) is locally 
smooth. However, the regions inside a closed leaf ip , where W is highly oscillatory, are of utmost 
importance for they neatly exhibit the nonpositive and thus nonclassical aspect of W which accounts 
for quantum interference and coherence phenomena. Thus, (0) provides an inadequate description 
of the propagation of Wigner functions at the semiclassical level. 

Most of these problems with the propagation of semiclassical Wigner functions were already clearly 
discussed by Heller M in 1976, even before the semiclassical approximations for W were properly 
developed 011]. The asymptotic expansion ([!]) in h was then partially resummed to obtain an 
improved propagation, though in a not very general context. From another approach, starting 
with the integral expression of the evolution equation, Marinov derived in 1991 a path integral 
representation for the propagator of Wigner functions, obtaining its semiclassical limit, which leads 
to (§). In this paper we develop the geometrical explanation of why the semiclassical limit of the 
propagation of Wigner functions cannot be derived from the semiclassical limit of their propagator. 

Thus, in §2 we present the integral equation for the Wigner evolution, defining its kernel, the 
Wigner propagator, whose stationary phase evaluation necessarily leads to classical propagation. 
Then, in §3 we recollect the constructions of semiclassical Wigner functions which will permit, in §4, 
a correct stationary phase evaluation of the full integral equation for their evolution. This leads, in 
§5, to a simple geometrical prescription for the semiclassical limit of the propagation of a Wigner 
function in terms of the classical flow of the tips of the chords centered on the arguments of the initial 
semiclassical Wigner function. We conclude in §6 with a discussion on the geometrical meaning of 
the difference between classical Liouville and semiclassical Wigner propagation. 

2. INTEGRAL EVOLUTION AND SEMICLASSICAL PROPAGATOR 

The starting point of our analysis is the expression for the product of operators in the center (or 
Weyl) representation. Thus, if A(x) and B(x) are the center representations of A and B , then the 
center representation for AB is given [§||10[ by the integral Moyal product: 



(4) [AB](x) = J dx'dx" A(x')B{x")exp(iA(x, x' , x") /H) , 

where A(x,x',x") = 2(x A x' + x' A x" + x" A x) is the symplectic area of the triangle with these 
midpoints. Iterating (|]) we obtain the integral equation for the evolution of Wigner functions as the 
Weyl transform of the quantum evolution equation p t = U- t poUt '■ 

(5) W t {x) = [ dx'dx"dx"'U. t {x')W Q {x")U t {x"')exp{2i(x Ax' + x" A x"')/K)8(x - x' + x" - x'") , 



where U t (x) is the Weyl propagator, i.e. the Weyl transform of the unitary evolution operator U t . 
The 5-function prescribes the 4 points in ]R 2n as vertices of a parallelogram. The phase is twice 
the area of this parallelogram (times Tr l ). But this is also the area of any element of a continuous 
family of quadrilaterals circumscribed to the given parallelogram. In other words, we can identify 
the symplectic area of the parallelogram with vertices at ( OC j Ob . OC . OC 1 clS half of the symplectic area 
of any quadrilateral whose sides are centered on these points. This brings the product rule for three 
operators in line with that for two operators and, indeed, the product of any number of operators 
will depend on the corresponding circumscribed polygon flOfl . Thus, denoting by A4 the area of a 
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quadrilateral as function of its midpoints, we have that A 4 (x, x', x", x'") = 2{x Ax' + i" Ax'") and this 
function is well defined only on the subset D 3 C (iR 2n ) 4 , isomorphic to (]R 2n ) 3 , determined by the 
^-function. Fixing one (say, the first, x) of these points, we obtain a subset D 2 C (IR 2n ) 3 , isomorphic 
to (IR 2n ) 2 . Denoting its induced measure by d 2 (x',x",x"') , we can rewrite (H) as 

(6) W t {x)= [ ^(x^x , ^x // 0^-t(x0^o(^O^t(^' / )e^(^A 4 (x,x',x ,, ,x ,,/ )//i) • 

Similarly, we can reparametrize the parallelogram by identifying: x" = x$ , x' = (x + Xq)/2 — fi , 
x'" = (x + x )/2 + /i , in which case we can use (|5]) to get an expression for the Wigner propagator 
or kernel L t (x ,x) , via its defining formula W t (x) = J dx L t (x , x)W (x ) , in the following form: 

(7) L t (x ,x) = J d/i U- t ((x + x )/2 - fj)U t ((x + x Q )/2 + fj,)exp(2i(fi A (x - x ))/H) . 

Marinov has derived an explicit path integral representation for the Wigner propagator L t (xo, x). 
This can also be achieved by introducing in ([/]) the path integral representation (TTJ| for the Weyl 



propagator U t (x) . For this latter, there exists a well-established semiclassical limit |]6[] IflOf 

(8) U t {x) « B lt (x)exp {iftT 1 [S lt — E lt t] (x)} , 

where S yt (x) is the symplectic area between the classical trajectory 7 4 (determined by the Weyl 
hamiltonian H) and the chord centered on x, with E lt = H(^ t ) being the energy of this trajectory, 
and where, again, £> 7t is a slow-varying real amplitude function. For sufficiently small times we can 
guarantee the existence of a single trajectory. Eventually there may be bifurcations, in which case 
(|) must be replaced by a sum of similar terms with appropriate Morse indices ||10|| . 



The naive semiclassical limit for the propagation of Wigner functions can be seen as a consequence 
of trying to get a semiclassical limit for the Wigner propagator ([?]) itself by stationary phase, either via 
its path integral representation or directly by using (||). In any case we are lead to a semiclassical 
Wigner propagator which is indeed classical, i.e. of the singular form: 

(9) Ct(x ,x) w 5(x - (a:)_ t ) , 

where x — ► (x)_ t is the inverse of the classical hamiltonian flow of H. Obviously, @ implies (|2|). 

Therefore, in order to obtain the correct semiclassical limit for the propagation of Wigner func- 
tions, one must apply stationary phase arguments to the full integral equation (f|f|) for the Wigner 
evolution, using the correct semiclassical expressions for the Wigner functions themselves. 

3. SEMICLASSICAL WIGNER FUNCTIONS 

For one degree of freedom, the simple form of the semiclassical Wigner function (|3]) depends on 
there existing only a single chord centered on each point x , besides the fulfilling of the semiclassical 
condition itself, i.e. the area Sf(x) being large in comparison to Planck's constant H. These conditions 
hold for points not too close to a convex leaf (curve) ip , outside the cusped triangular curve well in 
the interior of ij) known as the Wigner caustic |3J . Inside the caustic there are three chords centered 
on each point and the semiclassical Wigner function becomes a superposition of contributions of the 
same form as (Q), one for each chord and its associated area. Along the Wigner caustic itself, two 
chords coalesce and the amplitude in (|3|), defined as: 

(10) A^x) = (2cj/vt) (2nh\x + A x_\y 1/2 , 

blows up. Here, to is the classical frequency of motion along the curve ip , whereas x± are the phase 
space velocities at the tips x± of the corresponding chord centered on x. The Wigner caustic is hence 
the locus of the centers of the chords between points on ip with parallel or anti-parallel tangents. The 
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caustic condition thus leads to the vanishing of the skew-product x + A x_ and since this also happens 
when x converges onto if), this latter can be considered as a separate branch of the Wigner caustic, 
i.e. the breakdown of (|3|) also takes place on if) itself. Berry |J derived a uniform approximation 
based on the Airy function that is also oscillatory and asymptotically equivalent to (|3]) inside if) , 
rises to a smooth maximum near to this curve and decays exponentially outside. 



The same picture holds for integrable systems with L degrees of freedom [n}|. The novelty is 
that the caustic still has codimension 1, so the invariant L-torus arises as a higher singularity or 
catastrophe of the Wigner caustic. Another important feature is that multiple chords may be centered 
on points that lie arbitrarily close to the invariant torus. Still, sufficiently far within the energy-shell 
we retrieve the asymptotic form P"T|P|| 

(11) WO) « M x ) C0S i S j( x )/ h ~ njir/A} , 

3 

where again the actions Sj(x) are bounded by the j-th chord centered at x and any arc on the 
quantized torus (again a lagrangian leaf) between the chord tips. The amplitudes are now given by 

(12) A 3 (x) = (2/tt) [(2nh) L \det{I+,I.}\}- 1/2 , 

and rij is the signature of the matrix in ([12]). Here, the 2L action functions Ip are defined in terms 
of the L action functions / as lf{x) = I(x ± £j/2). Note that if L = 1 (|T2]) is identical with (|T0D. 



No immediate generalization of the semiclassical expression for pure-state wigner function is avail- 
able for chaotic or general nonintegrable systems. However, the superposition of pure Wigner func- 
tions in a classically narrow energy range e determines the spectral Wigner function W(x, E, e) which 
has the now familiar semiclassical limit || |10| 

(13) W(x, E,e)^J2 Aj(x)e-^ /h cos [Sj(x)/h - 7i ] , 

j 

where again Sj(x) is the symplectic area between a chord centered on x and a classical path between 
the chord tips. In this case, this is an arc of a classical trajectory along the energy shell , traversed 
in the time tj. For x close to the energy shell, the chords are all small and there is one arc traversed 
in a small time and a succession of longer arcs winding around the energy shell. Moving the center 
x leads to the crossing of many Wigner caustics where pairs of chords coalesce and disappear or 
vice versa. Along these caustics, one can again implement uniform approximations which coincide, 
off-caustics, with the sum (|H|) in which the amplitude for each chord is 

(14) Aj(x) = 2 L+1 {(27r7i)(£LE/(%)|cte£[l + M 3 }\y l/2 , 

where Mj is the matrix for the linearized classical map near the j-th trajectory arc on a special 
(2L — 2)-dimensional section that is centro-symmetric with respect to x. 

Thus we have a persistent overall picture: superposition of rapidly oscillating functions with wave 
vectors £j = —J[dSj/dx] , where J is the standard symplectic matrix in 1R 2L , and amplitudes 
depending only on the relation between the velocity vectors at the tips of each respective chord. 



4. WIGNER EVOLUTION: FULL STATIONARY PHASE 

Having recollected the general form of the semiclassical expressions for the various kinds of Wigner 
functions, we are now in a position to correctly approximate (0-|D by stationary phase. The cru- 
cial point here is how to geometrically interpret formulas (^-0) in the semiclassical approximation. 
As mentioned earlier, the phase of the exponential in the integrand, which is twice the area of the 
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parallelogram (times ft -1 ), corresponds to the area of any quadrilateral circumscribed to this par- 
allelogram. Furthermore, recall from the previous section that the argument x of the semiclassical 
Wigner function (|^ JTT|jr^) corresponds to the center of the chord whose tips lie on the energy shell 
or the lagrangian leaf ip corresponding to the quantum state ^ , so that the phase of its oscil- 
latory factor corresponds to the symplectic area between ip and the chord centered on x (minus 
7r/4). On the other hand, the semiclassical Weyl propagator (0) also has a phase corresponding to 
the symplectic area between the classical trajectory j t and the chord centered on x (minus E^t). 
Therefore, all the semiclassical terms of the integrand in (0-0) have oscillatory factors whose phases 
depend on the chords centered on their arguments. This suggests interpreting the phase A4 of the 
integrator in (0-0) as the area of the circumscribed quadrilateral that fits all the pertinent chords 
discussed above. To be totally consistent, such fitting must take into account all the correct ori- 
entations, so that one should first dismember the Wigner function as W = (W + + W~)/2 , where 
W ± (x) ~ A(x)exp{±i(S (x) / h — 7r/4)} , and propagate each member separately, then add them up 
for the total evolution of W. The situation for a single degree of freedom is illustrated in Figure 1. 



Fig. La Fig- Lb 




Stationary phase condition for the propagation of Wigner functions. 
Fig. l.a : Single degree of freedom. The various chords of ijj", jt and 7_j 
centered on x" , x, x'" and x' , respectively, match perfectly with the quadrilateral 
whose sides are centered on these points to yield, via eq.(16), the area between ip 
and the chord centered on x, the phase of the propagated Wigner function at x. 

Fig. l.b : The same situation but with all the orientations reversed. 

The really important issue is that these perfect matchings correspond precisely to the stationary 
phase condition, as we now show. Let us concentrate on the configuration in Fig. l.a, corresponding to 
the semiclassical evolution of Wq , determined by the leaf ip" , via formulas (0-0). The same analysis 
holds for the other case (Fig. l.b) as well as more dimensions. Inserting expressions (0) and (§) for 
Wo and U t in (0), we obtain: 

(15) W t + (x)=J ^(x^x^x , '0i3Xx0^ /, (x ,, )i3 /,, (x' / ')exp|^$ + (x,x',x /, ,x' ,/ )| , 
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$+(x, x', x", x'") = [S 7 _ t + E 7 _ t t]{x') + S r {x") - nh/A + [S 7t - E 7t t](x'") + A 4 (x, x' , x" , x'") . 

The stationary phase condition: <9<3> + /dx' = d& + /dx" = <9$ + jdx'" = implies that 

d[S^_ t + E 7 _ t t] dA 4 dS^ dA 4 d[S 7t - E 7i t] dA 4 
dx 1 dx' ' dx" dx" ' dx" 1 dx"' 



These equations mean |10| that the side of the quadrilateral A4 whose midpoint is x" coincides with 
the chord of S^n centered on x", with the same orientation as if>" . This is the chord from a to b in 
Fig(l.a). Similarly for the other cases, the side of A4 centered on x' is the chord connecting d to a 
and the side of A 4 centered on x'" is the chord connecting b to c. Therefore, the stationary phase 
condition coincides with the perfectly matching scenario. Furthermore, note that (E 7t — E 7 _ t )t is the 
area of the "curvilinear" quadrilateral which is formed by the paths if)" , 74 , if) -1 , 7_t , where if) is 
the immage of if)" under the hamiltonian flow for a time t. It follows by direct geometrical inspection 
on the perfectly matching configuration that 



(16) { S 7 _ t (x') + S r {x") + S 7t (x"') + A 4 (ar, x\ x", x'") - [E 7t - E 7 _ t )t = S^x) } 



match 



and thus $ttat( x ) = S A X ) ~ 71 h / 4 ■ And of course, $ stat = -$ttat ■ 

Therefore, we retrieve a new Wigner function W t whose phase is of the same general form given in 
(0,0,0) for Wt , in terms of the new action S^(x). Of course, this is just as we should expect within 
the general rules of semiclassical self-consistency in which all integrations are carried out within 
the stationary phase approximation. It should now be feasible to proceed with the full stationary 
phase evaluation of (0) for each of the different types of Wigner functions (|^JTT|JT3|) previously 



described, using the semiclassical expression H [10 for the amplitude B of the Weyl propagator 



Such a thorough verification of semiclassical self-consistency has already been carried out || for the 
pure state condition p 2 = p. It relies on such complicated geometrical constructions that we consider 
more appropriate, at this point, to just assume semiclassical self-consistency for the propagated 
amplitudes, as well. Accordingly, each semiclassically propagated Wigner function Wt retains its 
respective form (H,0,0), with its corresponding amplitude function (|IO|p2|,|T4|) . 



5. SEMICLASSICAL PROPAGATION: GEOMETRICAL PRESCRIPTION 

The beautiful consequence of this natural result is that it is sufficient to propagate the two tips 
of the chord of the original Wigner function Wo , centered on x" = x = (x + x$ )/2 , in order 
to obtain the value of the new Wigner function Wt at x = x t = (xf + xf )/2 . Then, using ( |T6|) 
we obtain the new phase S^x) = S^ t {x t ) from S^n{x") = S^ {xq). Similarly for the amplitudes, 
the actions S lt {x"') = {S 7t — E 7t t}(x" f ) and S 7 _ t (x') = {S 7 _ t + E 7 t t}(x') determine the matrices 
M + (x"') and M_(x') for the linearized classical maps from Xq = a;" = x'" to x'" = x + = xf and 
from %t = X- = x'_ to x' + = x"_ = Xq , respectively, by JTO] 



(17) M ± (y ± ) = [1 - J(d 2 S 7±t (y ± )/dyi)][l + J(d 2 S 7±t (y ± ) / 'dyl)]' 1 , 

where y + = x'", y_ = x'. Hence, the corresponding phase-space velocity vectors for the tips of the 
propagated chord are obtained from the original ones by 

(18) x+ = M + x[ ; x- = MZ x x"_ , 

which gives the new amplitudes A^(x) = A^ t (xt) from A^"(x") = A^ (xq) via their defining equa- 



tions ([T0|JT^JT^ ) . When the central action S 7t (y + ) or S 7 _ t (y^) goes through a caustic, one may rely 
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on the chord actions <S 7t (£+) or <S T _ t (£_) instead, for £± = ±(x^ — a: ), where the chord actions are 
the Legendre transform of the central actions: «S(£) = J-(y(£), £)> J~(y, £) = £ A y — S(y) , to get (10 



2 m-1 



(19) m ± (£ ± ) = -[i + j(a 2 <s 7±t (e ± )/^)][i - j(^<s 7±t (e ± )/^)] 

Note, however, that in this case one must be very careful in correctly counting the Maslov indices. 

Therefore, our analysis clearly leads to interpreting the correct point-to-point propagation of a 
semiclassical Wigner function as given by the simple geometrical picture: 

(20) W (x ) — > Wt(x t ) 

where x t is the midpoint of (x^, xf) , while (xq , Xq) stands for the tips of the original chord centered 
on xq . This prescription based on the tips-of-the-chord flow provides a precise semiclassical evalu- 
ation of W t (x t ) from Wo(x ). Specifically, QIBP is used to propagate the phase of W along the path 
Xq — > Xt ■ For short times, the smooth amplitude may be taken as approximately constant along 
this path, as illustrated below. For longer times the amplitude must also be properly propagated 



along this path via (OtHI)- in this way, fl20[ ) is straightforwardly used to obtain a quantitatively pre- 
cise semiclassical propagation of W "mostly everywhere". But, wherever Wt goes through an inner 
Wigner caustic, its analysis must be carefully completed by uniform approximations. Even then (p0|) 
provides most valuable information because ( p^ - |l9| ) can be used to determine these inner Wigner 
caustic points and ( |I6|) determine the contributing new phases which, from the knowledge of the 
corresponding caustics and uniform approximations for Wo and after a careful qualitative analysis 
of the flow ipQ — > ip t , can lead to the correct (Airy, Pearcey...0) functions representing Wt in these 
regions. The on-shell Wigner caustics are quite simpler to deal with because there x^ = x\ and ( p0|) 
coincides with (fj). Thus, we can say that (|20"D is quantitatively precise "mostly everywhere" and 
qualitatively precise everywhere, at the semiclassical level. 

As a simple partial illustration of the method, consider the hamiltonian h = (q 2 +p 2 ) 2 /4. Operator 
ordering is not relevant, modulo a constant, and H = h = r 4 /4 . The trajectories are circles around 
the origin, but 9 = r 2 so the classical flow of xq = (r ,9 = 0) is x t = {r t = r ,9 t = r^t). If xq is 
the center of the chord whose tips, lying on ipo , are Xq = (r_,a) , Xq = (r + , — then the new 
center x t = (r t ,9 t ) of (x^jxf) is given by f 2 = {r 2 _ + r^_)/4 + r^r+cos^ir 2 ^ — r 2 _)t — (a + (3))/ 2 , 
9t = r 2 _t + a + a' , with 2f t cos(a') = r_ + r + cos((r+ — r 2 _)t — (a + (3)) . Thus, not only 9 t ^ 9 t , but 
also r\ — r\ = r\ — r 2 , = r + r_{cos((r 2 + — r 2 J)t — (a + 0)) — cos(a + /3)}/2. 

Despite the oddity of this flow, Wt(x t ) is mostly everywhere obtained from Wo(xq) in a straight- 
forward manner from the knowledge of xf . Hence, modulo nontrivial Maslov changes, the phase 
difference 5S(x ,x t ) = S^ t (x t ) — S^ (x ) independs on ip or Wo and is (via (pT6|) ) given by: 

(21) 5S(x , x t ) = t{r% - ri)/4 + r + r_{sm{t{r 2 + - r 2 _)/2)cos(t(r 2 + - r 2 _)/2 - (a + (3))} . 

On the other hand, the new amplitude depends on the velocities of the tips of the chord deter- 
mined by ipo and thus depends on ipo and Wo itself in an intrinsic way. However, for short enough 
times such that x t is not too different from xt , we can approximate the slow-varying amplitude by 
At(xt) ~ At(xt) ~ ^4o(^o) , "using" (§) in a more justified way. Together with (^l|) this gives a first 
approximation for the semiclassical propagation of any Wigner function, mostly everywhere and for 
short times, under the r 4 hamiltonian. For longer times, we obtain At(x t ) from ^(^o) y i a (|10|) and 
(lilf), with M± determined by (\Tf\) from S± t = ±r|[r^t — 2sin(r±t)]/4 with \y±\ 2 = r 2 k cos 2 (r±t/2), 
and by (^) from the corresponding S±t . As for formula (0), again the precise error will depend on 
tpo and xq , but if xq is far from the leaf i/jq , we can approximate S^ t (x t ) — S^ (x ) oc t(r + — r_) 3 . 
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Thus, for any pair (ip , x ) such that r + and r_ differ significantly, (|2]) can be considerably wrong 
even for a short time propagation. 

6. CONCLUSION 

We have shown that the propagation of a semiclassical Wigner function is correctly determined by 
the classical flow of the tips of the chord whose center is the argument of the initial Wigner function, 
not by the classical flow of the argument itself. This reveals the irrelevance of the semiclassical 
limit for the Wigner propagator. Accordingly, in evaluating the integral equation for the evolution 
of a Wigner function, we found that satisfying the stationary phase condition was tantamount to 
matching four areas around an appropriate quadrilateral. However, by reducing the propagation 
to a single trajectory, traversed in positive and negative time, the side of the quadrilateral facing 
the chord of the Wigner function shrinks to a point and so it could only be matched by a zero 
length chord. Thus we found that the relevant trajectories for a precise semiclassical propagation 
are explicitly determined by the specific semiclassical Wigner function being propagated. 

Indeed, one way to correctly evolve a semiclassical Wigner function Wo determined by a leaf ipo in 
phase space is to classically evolve ipo via the hamiltonian flow: tpo — > ip t and then evaluate the new 
semiclassical Wigner function W t at each point from the knowledge of ip t via (|3|, |T0| - [r4|) , as shown 



in 0. Thus, to obtain the value of the new semiclassical Wigner function Wt at x one needs to 
know, at the very least, the corresponding chord of ipt centered on x. In fact, our analysis shows 
that the knowledge of the flow of each pair of points in ipo is sufficient to determine the evolved 
semiclassical Wigner function Wt at all points. Though these two prescriptions for the evolution of 
W are equivalent, in practice the semiclassical point-to-point propagation of W is often simpler to 
determine, both qualitatively and quantitatively, by ( ]TE| - |r9| ) via tips-of-the-chord flow (fZH|). 

On the other hand, a mere knowledge of the hamiltonian flow of an argument of Wo is not enough 
to reconstruct the new chord and hence the new value of Wt- The exceptions are either linear flow 
or when the two tips of the corresponding chord coalesce into its center, in other words, when the 
argument of Wo lies on the leaf ipo . For these points the total chord flow coincides with the center 
flow, thus (@) is verified in the regions of high amplitude near the leaf ipo, but not otherwise for 
nonlinear flows. It is easy to see why: the phase of Wq(xo) corresponds to the symplectic area 
between ipo and the chord centered on xq. Under a nonlinear hamiltonian flow, the symplectic area 
of this closed circuit is preserved, however the immage of the chord is no longer a straight segment 
and thus the symplectic area between ip t and a new chord centered on x t will generally differ from 
the corresponding area evaluated at t — 0. This difference can be large, even for short times, if 
the tips of the initial chord are sufficiently far apart, i.e. when the argument of the initial Wigner 
function is sufficiently far from the initial classical leaf. Therefore, generally (fj) is simply wrong, 
semi classically. In sharp contrast, the prescription based on the tips-of-the-chord flow (|20| ) provides a 
precise semiclassical evaluation of Wt(x t ) from Wq(xq). Only when x t = Xt do (0) and fl20|) coincide. 
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